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Abstract 

We investigate the inclusive production of prompt J /if) mesons in polarized 
hadron-hadron, photon-hadron, and photon-photon collisions in the factorization 
formalism of nonrelativistic quantum chromodynamics (NRQCD) providing all con- 
tributing partonic cross sections in analytic form. In the case of photoproduction, we 
also include the resolved-photon contributions. We present numerical results appro- 
priate for BNL RHIC-Spin, the approved SLAC fixed-target experiment E161, and 
the e + e~ and 77 modes of TESLA. Specifically, we assess the feasibility to access the 
spin-dependent parton distributions in the polarized proton and photon. We also 
point out that preliminary data on J/ip inclusive production taken by the PHENIX 
Collaboration in unpolarized proton-proton collisions at RHIC tends to favor the 
NRQCD factorization hypothesis, while it significantly overshoots the theoretical 
prediction of the color-singlet model at large values of transverse momentum. 

PACS numbers: 12.38.Bx, 13.60.Le, 13.85.Ni, 14.40. Gx 



1 Introduction 



Since its discovery in 1974, the J/ip meson has provided a useful laboratory for quantita- 
tive tests of quantum chromodynamics (QCD) and, in particular, of the interplay of per- 
turbative and nonperturbative phenomena. The factorization formalism of nonrelativistic 
QCD (NRQCD) 1J provides a rigorous theoretical framework for the description of heavy- 
quarkonium production and decay. This formalism implies a separation of short- distance 
coefficients, which can be calculated perturbatively as expansions in the strong-coupling 
constant a s , from long-distance matrix elements (MEs), which must be extracted from 
experiment. The relative importance of the latter can be estimated by means of velocity 
scaling rules, i.e., the MEs are predicted to scale with a definite power of the heavy-quark 
(Q) velocity v in the limit v <C 1. In this way, the theoretical predictions are organized 
as double expansions in a s and v. A crucial feature of this formalism is that it takes into 
account the complete structure of the QQ Fock space, which is spanned by the states 
n = 2S + l l}j> with definite spin S, orbital angular momentum L, total angular momentum 
J, and color multiplicity c = 1,8. In particular, this formalism predicts the existence of 
color-octet (CO) processes in nature. This means that QQ pairs are produced at short 
distances in CO states and subsequently evolve into physical, color-singlet (CS) quarkonia 
by the nonperturbative emission of soft gluons. In the limit v — > 0, the traditional CS 
model (CSM) |2|3| is recovered. The greatest triumph of this formalism was that it was 
able to correctly describe [3] the cross section of inclusive charmonium hadroproduction 
measured in pp collisions at the Fermilab Tevatron [5] , which had turned out to be more 
than one order of magnitude in excess of the theoretical prediction based on the CSM. 

Apart from this phenomenological drawback, the CSM also suffers from severe concep- 
tual problems indicating that it is incomplete. These include the presence of logarithmic 
infrared divergences in the 0(a s ) corrections to P-wave decays to light hadrons and in 
the relativistic corrections to S"-wave annihilation [0] , and the lack of a general argument 
for its validity in higher orders of perturbation theory. While the fc^-factorization j7j and 
hard-comover-scattering |Hj approaches manage to bring the CSM prediction much closer 
to the Tevatron data, they do not cure the conceptual defects of the CSM. The color 
evaporation model which is intuitive and useful for qualitative studies, also leads to 
a significantly better description of the Tevatron data, but it is not meant to represent 
a rigorous framework for perturbation theory. In this sense, a coequal alternative to the 
NRQCD factorization formalism is presently not available. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and uN deep-inelastic 
scattering (DIS), e + e~ annihilation, 77 collisions, and 6-hadron decays may be found in 
the literature; see Ref. ^Uj and references cited therein. Furthermore, the polarization 
of charmonium, which also provides a sensitive probe of CO processes, was investigated 
|ll|12j . Until very recently, none of these studies was able to prove or disprove the 
NRQCD factorization hypothesis. However, HI data of ep — > eJ/ip + X in DIS at HERA 
[T3] and DELPHI data of 77 —>■ J/ip + X at LEP2 ^3] provide first independent evidence 
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for it [THUKj . 

Before the advent of the NRQCD factorization formalism, heavy- quarkonium produc- 
tion in experiments with polarized proton or photon beams was believed to provide reliable 
information on the spin-dependent gluon distribution functions of the polarized proton or 
photon. At present, however, we are faced with the potential problem that NRQCD pre- 
dictions at lowest order (LO) come with a considerable normalization uncertainty. This 
is due to the additional expansion in v, whose convergence property has yet to be tested, 
and the introduction of CO MEs as additional input parameters. In fact, the relevant 
CO MEs have only been determined through LO fits to experimental data and thus come 
with appreciable theoretical uncertainties. Thus, it must be clarified if heavy- quarkonium 
production with polarized proton or photon beams remains to be a useful probe of the po- 
larized gluon distribution functions. It is the purpose of this paper to answer this question. 
Specifically, we consider inclusive J /if) production in polarized pp, jp, and 77 collisions, 
appropriate for RHIC-Spin, experiment E161, and the DESY TeV-Energy Superconduct- 
ing Linear Accelerator (TESLA) operated in the e + e~ and 77 modes, respectively. In the 
RHIC-Spin mode of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory (BNL), proton beams with strong longitudinal polarization, of approximately 
70%, collide with center-of-mass (cm.) energy up to yS = 500 GeV and luminosity 
C = 2 x 10 32 cm~ 2 s -1 |17j . In the approved experiment E161 at the Stanford Linear 
Accelerator Center (SLAC), circularly polarized photons with energies between 35 and 
48 GeV will collide on a fixed target made of longitudinally polarized deuterium [T81. 
At TESLA, which is being designed and planned at the German Electron Synchrotron 
(DESY) in Hamburg, longitudinal electron and positron polarizations of up to 80% and 
60%, respectively, can be achieved. In the e + e~ mode [TH], photons are unavoidably 
generated by hard initial-state bremsstrahlung and by beamstrahlung, the synchrotron 
radiation emitted by one of the colliding bunches in the field of the opposite bunch. In 
both cases, polarization is transferred from the electrons and positrons to the photons in 
a calculable way. TESLA can also be converted into a photon collider via Compton back- 
scattering of high-energetic laser light off the electron beams [201 ■ With 100% circular 
laser polarization and up to 80% longitudinal electron polarization, the scattered photons 
can be arranged to have a strong circular polarization, especially at the upper and lower 
ends of their energy spectrum. 

The photons can interact either directly with the quarks participating in the hard- 
scattering process (direct photoproduction) or via their quark and gluon content (resolved 
photoproduction). Thus, the process 77 — > J/ip + X, where X collectively denotes all 
unobserved particles produced along with the J/ip meson, receives contributions from 
direct, single-resolved, and double- resolved channels. All three contributions are formally 
of the same order in the perturbative expansion and must be included. This may be 
understood by observing that the parton distribution functions (PDFs) of the photon 
have a leading behavior proportional to a In (ftj / Aq CD ) oc a/a s , where a is Sommerfeld's 
fine-structure constant, /i/ is the factorization scale, and Aqcd is the asymptotic scale 
parameter of QCD. Similarly, the reaction ^yN — > J/ip+X, where N is a nucleon, proceeds 
via direct and resolved channels, which must all be taken into account. 
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The J / ip mesons can be produced directly; or via radiative or hadronic decays of heav- 
ier charmonia, such as XcJ an d if}' mesons; or via weak decays of b hadrons. The respective 
decay branching fractions are B(xco — ► J/ip + 7) = (1-02 ± 0.17)%, B(xci J/4> + 7) = 
(31.6 ± 3.2)%, B( X c2 -> J/ij) + 7) = (18-7 ± 2.0)%, B(if/ -> J/ip + X) = (55.7 ± 2.6)%, 
and B{B -> J/ip + X) = (1.15 ± 0.06)% J23- At RHIC-Spin and TESLA, the b hadrons 
can be detected by looking for displaced decay vertices with dedicated vertex detectors. 
At RHIC-Spin, only about 1% of all J/ip mesons are expected to originate from 6-hadron 
decays [22]- I n experiment E161, the cm. energy is not sufficient to produce bb pairs. In 
our analysis, we will not consider J/ip mesons from 6-hadron decays. The cross sections 
of the four residual indirect production channels may be approximated by multiplying the 
direct-production cross sections of the respective intermediate charmonia with their decay 
branching fractions to J/ip mesons. It has become customary to collectively denote the 
J /if) mesons produced directly or via the feed-down from heavier charmonia as prompt. 

This paper is organized as follows. In Sec. |2l we list our formulas for the contributing 
partonic cross sections and compare them with the available literature results. Lengthy 
expressions are relegated to Appendix. In Sec. El we present our numerical results and 
discuss their phenomenological implications for the three experiments under consideration. 
Our conclusions are summarized in Sec. 0] 



2 Analytical results 

Let us consider the generic inclusive process AB —> H + X , where A and B are the 
beam particles (hadrons or photons) with definite helicities £a and respectively, and 
H is the charmonium state, and let a^ A ^ B denote its cross section. Then, the double 
longitudinal-spin asymmetry is defined as 



Act . . 

All = , (1) 

a 



where 
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Aa = - (-1)* A * B o*** 3 (2) 

are the unpolarized and polarized cross sections, respectively. An alternative definition 
of All, which allows one to study the dependences on the transverse momentum and 
rapidity y of H, uses the differential cross sections dAa and da. 

Invoking the factorization theorems of the QCD parton model and NRQCD, the dif- 
ferential polarized cross section can be written as 

dAa(AB ^H + X) = E / dx a Af a/A (x a , (i f ) f dx b Af b/B (x b , /i f ) ^(0 H [n\) 

a,b,d J J n 

x dAa(ab — > cc[n]d), (3) 
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where it is summed over 0,6 = 7, 9i <2S Q an d d = g,q, q, x a is the fraction of longitudinal 
momentum that a receives from A, fif is the factorization scale, (O [n]) are the MEs of 
the if meson, and 



dAa(ab — > cc[n]d) = — do + (ab — > cc[n](i) — da~(ab — ► ce[n]d) 



(4) 



Here, f^ A (x a , fif) are the PDFs of A for the case that the spin of a is (anti)parallel to that 



of A (£a£o, = an d da ± (ab — > cc[n]d) refers to the case of £ a £& = ±1. In fact, due to 
parity conservation in QCD, da ± (ab — > cc[n]d) only depends on the product as can be 



seen from the Appendix. With the definition A/ 7 / 7 (x 7 , p,f) = 5(1- 



7/' 



Eq. © 



accommo- 



dates the direct, single-resolved, and double-resolved channels. In the case of TESLA, ad- 
ditional convolutions with appropriate photon flux functions Af^/ e (x) = f~u{x) — /, 



lj/e 



x) 



where the superscript ± refers to the case of £ e £ 7 = ±1, have to be implemented in 
Eq. (|5]L invoking the Weizsacker- Williams approximation (WWA) [23] for electromag- 
netic bremsstrahlung and its analogues for beamstrahlung and Compton back-scattering. 
The corresponding formulas for the unpolarized case are obtained by omitting every- 
where in Eq. (JHJ) the A symbols. The kinematical relations needed to calculate pt and y 
distributions from Eq. (J3J may be found, e.g., in Ref. 

We work in the fixed-flavor-number scheme, i.e., we have rif = 3 active quark flavors 
q = u, d, s in the proton and resolved photon. As required by parton-model kinematics, 
we treat the quarks q as massless. The charm quark c and antiquark c only appear in 
the final state. The cc Fock states contributing at LO in v are n = ^(p, ^jf , 3 S[ 8 \ 

lp(8) 



if H 



„ . „ _ 3o(l) 10(B) 30(8) 

and n = 3 P?\ 3 si 8) if H = Xc j 
relations 



3p 



(8) 
J 



if if = 

where J 



JM V, ^(35), ■ ■ n = x Pr, %*> if H = h c ; 
= 0,1,2. Their MEs satisfy the multiplicity 



Qip(nS) 
QXcJ 
QXcj 



3p(8) 
3p(l) 
'35.(8) 



(2J + 1) (C^™ 5 ) 



3p(8) 



(2J + 1) (O Xc0 



(2J + l)(O Xc0 ^1 



3p(l) 
3o(8)' 



(5) 



which follow to LO in v from heavy-quark spin symmetry. In our numerical analysis, 
we only include the J/ip, Xcj, and ip' mesons. For completeness and future use, we also 
list formulas for all the other experimentally established charmonia. Our results readily 
carry over to all known bottomonia as well. In order for the J /if) meson to have finite 
transverse momentum px, we allow for an additional parton (quark or gluon) in the final 
state. We do not include the J/x/j + j final state because prompt photons can be identified 
experimentally, so that such events can be eliminated from the data sample. 
We are thus led to consider the following partonic subprocesses 
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cc[^]g, 
cc[^]g, 



(6) 
(7) 
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ccK 
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(»: 



(9) 
(10) 

(ii) 

(12) 
(13) 
(14) 
(15) 



where q = l S , 3 Si, 1 Pi, 3 Pj with J = 0,1,2. The processes 77 — > cc[^}g are forbidden 
by color conservation. Similarly, the processes jq — > cc[^]q are prohibited because the 
c-quark line is connected with the g-quark line by one gluon, which transmits color to 
the cc pair. By angular-momentum conservation, processes (0) and (J7|) with ^ = 1 So, 3 Pj, 
process © with q = x Pi, and processes (|T2"j) and (fTl^) with ^ = 3 S±, X P\ have zero cross 
sections. We calculated the differential cross sections da/dt of all remaining partonic 
subprocesses for arbitrary helicities £ a ,6> = ±1 of the incoming particles a,b = l,g,q,q. 
Our formulas are listed in the Appendix. 

In the literature, inclusive charmonium photoproduction with polarized beams was 
studied both in the CSM |25l2(i| and in NRQCD |27l28l29j . In the case of hadroproduction, 
CSM and NRQCD studies may be found in Refs. |3f26p30p3"T| and Ref. [3*2] . respectively. 
We are not aware of any previous analyses for collisions of polarized photon beams. Some 
of these papers contain analytic results |3p25p2(j|27p28|30j . which we can compare with ours. 
This is necessary because the literature results are in some cases mutually conflicting. As 
for direct J / ip photoproduction, our result for the CS process (J7J) with ^ = 3 S% agrees with 
Refs. |2fil27j . while it disa grees with Ref. [2H] ; our results for the CO processes (JBI) and © 
with <; = ^o, 3 Si, 3 Pj agree with Ref. (2Z] 1 and partly with Ref. [2H]- 2 As for charmonium 
hadroproduction in the CSM, our results for processes (fiT)|) agree with Refs. [3 30J and 
those for processes (fl"2"|) and (fHI|) with q = 3 Pj agree with Ref. [SU]. To our knowledge, 
our residual formulas are not available in the literature. 

3 Numerical results 

We are now in a position to present our numerical results. We first describe our theoretical 
input and the kinematic conditions. We use m c = (1.5 ± 0.1) GeV, a = 1/137.036, and 
the LO formula for a^ f \fj>) |2I] with rif = 3 active quark flavors and asymptotic scale 
parameter A^ 3 ) = 204 MeV. Here, m c denotes the charm-quark mass. 

1 We do not agree with Eqs. (A6) and (A7) of Ref. [22] and find an overall minus sign relative to 
Eq. (All). However, we agree with Eqs. (A12) and (A13), which refer to the sum over the three 3 Pj 8 ^ 
states. 

2 We agree with Eq. (A3) of Ref. [2H if we divide the right-hand side of this equation by s 2 . Equa- 
tions (A4)-(A6) are corrupted, contain undefined symbols, and are thus unsuitable for comparisons. 
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As for the polarized proton PDFs, we use the LO sets by Gliick, Reya, Stratmann, and 
Vogelsang (GRStV) [33] and those by Gehrmann and Stirling (GS) jMj. Their reference 
sets of unpolarized proton PDFs are both by Gliick, Reya, and Vogt, namely the LO sets 
of Refs. [351 (GRV98) and [36j (GRV95), respectively. There are two GRStV sets, one 
of which corresponds to the common standard (STD) scenario of polarized PDFs with a 
flavor-symmetric light-sea (antiquark) distribution Afg/ P , while the other one refers to a 
completely 577(3) /-broken valence (VAL) scenario with totally flavor-asymmetric light-sea 
distributions (Afe/p ^ Afj, p ^ Afs/ p ). The latter are modelled with the help of a Pauli- 
blocking ansatz at the low radiative or dynamical input scale [i^o = 0-26 GeV 2 , which 
complies with predictions of the chiral quark-soliton model and expectations based on the 
statistical parton model as well as with the corresponding flavor-broken unpolarized sea 
if lip > fu/p)- There are three GS sets (A,B,C), corresponding to three different, equally 
possible scenarios for the polarized gluon distribution, which is only loosely constrained 
by current experimental data. Leaving aside nuclear corrections and appealing to strong- 
isospin symmetry, the effective nucleon (JV) PDFs may be approximated by 



f a/N (x,n f ) = — \Zf a/j ,(x,fif) + (A- Z)f b/P (x,fif) 



(16) 



and similarly for Af a /x{x, ///), where Z and A are the respective numbers of protons 
and nucleons in the nucleus and (a, b) = (it, d), (u,d), (d,u), (d,u), (s,s), (s,s), (g,g). In 
the case of deuterium D, we have Z — 1 and A — 2. We present the cross sections per 
nucleon, rather than per nucleus. 

Up-to-date LO sets of polarized photon PDFs were presented by Gliick, Reya, and Sieg 
(GRSi) [37] . Their unpolarized counterpart is the LO set by Gliick, Reya, and Schienbein 
(GRSc) 38 . There are two GRSi sets, which differ in the assumed boundary condition 
for the polarized gluon distribution at the dynamical input scale and are characterized 
as maximally (MAX) or minimally (MIN) saturated. Our default sets are taken to be 
GRStV-STD and GRSi-MAX as well as their unpolarized counterparts GRV98 and GRSc, 
respectively. The other sets are employed to assess the sensitivity of the considered 
experiments to the details of the polarized proton and photon PDFs given the theoretical 
uncertainties from other sources, especially those introduced by NRQCD. 
We choose the renormalization and factorization scales to be //j 



ii\jkm 2 c +p?t, with 

i — r,f, respectively, and independently vary the scale parameters £ r and £/ between 
1/2 and 2 about the default value 1. As for the J/ip, Xcj, and ip' MEs, we adopt the 
set determined in Ref. [12] using the LO proton PDFs by Martin, Roberts, Stirling, and 



Thorne [H]. Specifically, {0^ nS ^ 3 S[ 



with n - 

from the measured partial decay widths of ip(nS) 
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while ( 0^ nS ^ 



15.(8) 



Qip(nS) 3C7j8) \ /q^{uS) 



3p(8) 



1, 2 and (O x * 
• l + l~ and Xc2 ~ 
and ( O x -° 



3 P, 



(i) 


77 
35(8) 



were extracted 
2"Tj . respectively, 
were fitted to 



the px distributions oiip(nS) and Xcj inclusive hadroproduction [5] and the cross-section 



ratio a x Ja 



Xcl 



measured at the Tevatron. The fit results for (C^( n5 ) 



15.(8) 
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lQil>{nS) 



3 P f 



(8) 



are strongly correlated, so that the linear combinations 



j^ip(nS) _ /Qif>(nS) 



lo(8) 
°0 



_| / Qi>(nS) 



mi 



3 P f 



(8) 



(17) 



with suitable values of r, are quoted. Unfortunately, the hadronic cross sections under 



(8) 



consideration here are sensitive to different linear combinations of (o^^^ 



3 P, 



write ( 0^ nS) 



and 

than those appearing in Eq. (fTTj) . In want of more specific information, we 

15.(8) 



= K^ nS) Mf^ and (O^ [ 3 P {8) ]) = (1 - K f(nS) ) {m 2 Jr)Mf nS ) 
and independently vary kj/^ and k^i between and 1 around the default value 1/2. Our 
choice of A^ 3 ) coincides with the value employed in Refs. |33|35|37|38|39| . and it is very 
close to the value 232 MeV of Refs. pi3T)j . 

In order to estimate the theoretical uncertainties in our predictions, we vary the un- 
physical parameters £ r , £/, and k$> as indicated above and take into account the 

experimental errors on m c , the decay branching fractions, and the MEs. We then com- 
bine the individual shifts in quadrature, allowing for the upper and lower half-errors to 
be different. 

We now discuss the photon flux functions that enter our predictions for photon-photon 
scattering at TESLA. The energy spectrum of the bremsstrahlung photons is well de- 
scribed in the WWA. The formulas for the unpolarized and polarized cases, may be found 
in Eq. (20) of Ref. |H] and Eq. (14) of Ref. 02], respectively, where the maximum photon 
virtuality Q 2 is given in Eq. (26) of Ref. |H]. We assume that the scattered electrons and 
positrons will be antitagged, as was usually the case at LEP2, and take the maximum 
scattering angle to be max = 25 mrad jUJ]. The energy spectrum of the beamstrahlung 
photons is approximately described by Eq. (2.14) of Ref. jH] in the unpolarized case and 
by Eq. (7.10) of Ref. [IS] (see also Ref. jl^]) in the polarized one. It is controlled by the 
effective beamstrahlung parameter T, which is given by Eq. (2.10) of Ref. jH]. Inserting 
the relevant TESLA parameters for the = 500 GeV baseline design specified in Ta- 
ble 1.3.1 of Ref. ^H] in that formula, we obtain T = 0.053. In the case of the e + e~ mode 
of TESLA, we coherently superimpose the WWA and beamstrahlung spectra. Finally, 
in the case of the 77 mode of TESLA, the energy spectrum of the back-scattered laser 
photons is given by Eq. (6a) of Ref. [UJ in the unpolarized case and by Eq. (4) of Ref. 
in the polarized one. It depends on the parameter k = s n /m 2 — 1, where 
energy of the charged lepton and the laser photon, and it extends up to x max = k/(k+ 1), 
where x is the energy of the back-scattered photons in units of y/S/2. The optimal value 
of k is k = 2 (l + a/2^ ps 4.83 jUJ, which we adopt; for larger values of k, e + e~ pairs 
would be created in the collisions of laser and back-scattered photons. 

Our main numerical results are presented in Figs. ^ El El and EJ which refer to pp —>■ 
J/ip+X at RHIC-Spin, to 7-D — > J/ip+X in Experiment E161, and to e + e~ — > e + e~ J/ip+ 
X at TESLA in the e + e~ and 77 modes, respectively. In each figure, the NRQCD 
and CSM predictions for the unpolarized cross sections da/dpT and da/dy or da/dz are 
displayed in the first panel, while those for the double longitudinal-spin asymmetry All 
are shown in the second and third panels, respectively. For the collider experiments, we 



s e7 is the cm. 
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consider da/dy, while for the fixed-target scattering process jD — > J/tp + X, we consider 
da/dz, where z = {pj/^-Pd) /{v~i'Pd) = \J m2 j + Pt cosh y/E 7 is the inelasticity variable. 
In Fig. |2l which refers to the latter case, the first panel also contains NRQCD and CSM 
predictions for the resolved-photon contributions by themselves. In the left and right 
columns of each figure, the px and y or z dependences, respectively, are studied. As 
for RHIC-Spin, E161, and TESLA, the pt distributions are integrated over the intervals 
\y\ < 2.4, 0.3 < z < 0.8, and \y\ < 2, while the y or z distributions are integrated over 
all kinematically allowed values of px in excess of 15, 1.5, and 10 GeV, respectively. In 
Figs. ^ HU and 01 which refer to collider experiments, the y distributions are symmetric 
about y = 0. The shaded bands indicate the theoretical uncertainties in the NRQCD and 
CSM default predictions, excluding the freedom in the choice of the PDFs. In the case of 
All, these uncertainties are compared with the spread due to a variation of the polarized 
PDFs. We assume the ideal case of 100% beam polarization, except for TESLA, where we 
take P{e + ) = 60% and P(e~) = 80%. As for RHIC-Spin and E161, realistic polarization is 
straightforwardly accounted for by scaling All with [P(p)] 2 and P(7)P(P), respectively. 

Considering the unpolarized cross sections, we observe in all cases that the NRQCD 
predictions dramatically exceed the CSM ones, by typically one order of magnitude at 
small values of pt- This excess tends to amplify as px increases, and it is striking in the y 
distributions, thanks to the relatively large minimum-p^ cuts chosen. From the first panel 
of Fig. El we observe that the resolved-photon contribution to 7P —>■ J ftp + X is greatly 
suppressed in the entire kinematical range considered, both in NRQCD and the CSM. For 
the same reason and because the polarized photon PDFs are poorly known at present, 
we neglect this contribution in the evaluation of All presented in the second and third 
panels of Fig. 121 By contrast, the NRQCD and CSM predictions for TESLA, shown in 
Figs. 01 and El are greatly dominated by the single- resolved contributions, as was already 
observed in Ref. [4*9*] . We conclude that, in all these experiments, the normalization of the 
unpolarized cross section is a distinctive discriminator between NRQCD and the CSM. 

In fact, data from the PHENIX Collaboration [22] at RHIC, with v 7 ^ = 200 GeV, 
tend to favor the NRQCD prediction compared to the CSM one. This is demonstrated in 
Fig. 03 where the differential cross section d 3 a/dyd 2 px is analyzed for 1.2 < y < 2.2 as a 
function of px- Since the J/ip mesons are tagged through their decays to pairs, the 

factor B(J/ip — > = (5.88 ± 0.10)% [21] is included in the theoretical predictions. 

We observe that, for p T > 2 GeV, the data is nicely described by the NRQCD prediction, 
while they significantly overshoot the CSM ones. A similar observation was made in 
Ref. [50J. On the other hand, such a comparison has to be taken with a grain of salt in 
the bin 1 GeV < px < 2 GeV because, at LO, the NRQCD prediction and the P-wave 
contribution to the CSM one suffer from infrared and collinear singularities at px = 0, 
which still feed into that bin as an artificial enhancement. 

Looking at the second and third panels of Figs, d El and|H we observe that different 
choices of polarized proton and photon PDFs yield discriminative NRQCD and CSM pre- 
dictions, as far as the colliding-beam experiments are concerned. In fact, the differences 
are large against the combined theoretical uncertainties from all other sources, especially 
at large values of px- This means that sufficiently precise measurements of All in these 
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experiments are bound to increase our knowledge on the spin-dependent parton structure 
of the polarized proton and photon, regardless of the presently still somewhat unsatisfac- 
tory status of the residual theoretical uncertainties. As for the pt distribution of All at 
RHIC-Spin, the NRQCD and CSM predictions incidentally almost coincide, so that the 
polarized proton PDFs can be explored in a model-independent fashion. On the other 
hand, the NRQCD and CSM predictions for the y distribution of All at RHIC-Spin 
exhibit strikingly different shapes in the forward and backward directions. 

As for the fixed-target experiment E161, the situation is somewhat less favorable 
because the theoretical uncertainties not related to the polarized proton PDFs are larger, 
owing to the relatively low photon-nucleon cm. energy v^S* = m^(2E^ + mj^) « 9.2 GeV. 
This is especially the case for the NRQCD prediction shown in the second panel of Fig. |21 
Nevertheless, with enough experimental statistics, it should be possible to discriminate 
between the GRStV and GS sets of polarized proton PDFs. In the CSM, the resolving 
power of All with respect to the spin-dependent parton structure of the polarized proton 
used to be much better, as is evident from the third panel of Fig. El In this sense, 
the introduction of NRQCD, which was necessary to overcome phenomenological and 
conceptual problems of the CSM, led to some aggravation. 

4 Summary 

Using the NRQCD factorization formalism at LO, we studied the inclusive production 
of prompt J/ip mesons in polarized hadron-hadron, photon-hadron, and photon-photon 
collisions at RHIC-Spin, experiment E161, and TESLA, respectively, with regard to the 
potential of these facilities to resolve the spin-dependent parton structure of the polarized 
proton and photon. For future use by other authors, we provided all contributing partonic 
cross sections in analytic form. 

Our main message is that the spread in the theoretical predictions encountered by 
trying out in turn various up-to-date sets of polarized proton and photon PDFs in general 
considerably exceeds the combined theoretical uncertainties from other sources, which we 
estimated rather conservatively. Therefore, these experiments have the power to improve 
our knowledge of the spin structure of the proton and photon already at present. This re- 
solving power can be increased in the future by including next-to-leading-order corrections 
in a s and v. 

As a by-product of our analysis, we found that preliminary PHENIX data of unpo- 
larized hadropro duct ion of J/ijj mesons at RHIC tends to favor NRQCD as compared to 
the CSM. This is in line with previous findings in connection with hadroproduction at 
the Tevatron 0j, electroproduction at HERA [T3], and photoproduction at LEP2 [To*] . 

Acknowledgements 

We thank Peter Bosted and Gerry Bunce for useful communications ^H] and Christoph 
Sieg for providing us with a parameterization of the GRSi polarized photon PDFs [37 . 
This work was supported in part by the Deutsche Forschungsgemeinschaft through Grants 



10 



No. KL 1266/1-3 and KN 365/1-1, by the Bundesministerium fur Bildung und Forschung 
through Grant No. 05 HT1GUA/4, and by Sun Microsystems through Academic Equip- 
ment Grant No. EDUD-7832-000332-GER. 



A Partonic cross sections 

Here, we list the differential cross sections da/dt of processes ( j H|)-()15 |l for arbitrary he- 
licities £ a ,£fc of the incoming particles. The Mandelstam variables s, t, and u are defined 
in the usual way, and M = 2m c , where m c is the charm-quark mass. By four-momentum 
conservation, we have s + t + u = M 2 . Furthermore, e = \/A~kq. and g s = \J Aixct s are the 
electromagnetic and strong gauge couplings, and Q q is the fractional electric charge of 
quark q. Our results read: 



da f 
- (77 - cc 



35.(8) 



eYs (O 



35.(8) 



37TS 2 (s + tf(s + uf(t + uf 



X 



t 2 u 2 + stu (t + u) + s 2 (t 2 + tu + m 2 ) + £ a £ t 



— tu — stu (t + u) + s tu 



(A.l) 



da 
dt 



77 — > cc 



lp(8) 



x I tuit + uf (t 2 + tu + u 2 ) + s[f + Atu + lltV + HtV + lltV + Atu b + u 
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32e 4 g 2 (0 )Q. 



3irMs 2 (s + tf(s + uf(t + uf 



4„ 2 



i2„ 4 



■ 4 | 

2 



+ s 2 (t + u) (At 4 + 7t 3 u + 9t 2 u 2 + 7tu 3 + Au 

+ s 3 (jt 4 + Ut 3 u + 16tV + Utu 3 + 7m 4 ) + s 4 (t + u) (jt 2 + Atu + 7u' 
+ As 5 (t 2 + tu + m 2 ) + s 6 (t + u) + fV(t + uf + 2stw(t 2 + tu + u 2 ' 
+ s 2 {t + u) (t 4 + At 3 u + 3tV + Atu 3 + u 4 ) + 2s 3 (2t 4 + 5t 3 u + 5t 2 u 2 + 5tu 3 + 2u 

+ 2s 4 (t + u) (3t 2 + 2tu + 3u 2 ) + 4s 5 (t 2 + tu + u 2 ) + s 6 (t 



(A.2) 



da_ 
dt 



(l9 



cc 



3 sP 



Ac 2 g 4 s (0 fei 1} l)g 2 M 



97rs 2 (s + tf(s + uf(t + uf 



X 



t 2 U 2 + stM (t + U) + S 2 (t 2 +tu + M 2 ) + £ a £fc 



— t u — stu (t + u) + s tu 



(A.3) 



cla / 

dt (™ C ' 



ip(D 



16e 2 5 4 ((9 fe^DQ 2 



x 



" 9nMs 2 (s + tf(s + uf(t + uf 
\tu(t + uf (t 2 + tu + u 2 ) + s (t 6 + At b u + lltV + 14tV + lltV + Atu 5 + u 6 ) 
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+ s 2 (t + u) (At 4 + 7t 3 u + 9t 2 u 2 + 7tu 3 + 4m 4 ) 
+ s 3 (7t 4 + Ut 3 u + 16tV + Utu 3 + 7m 4 ) 

+ s 4 (t + u) (7t 2 + 4tw + 7m 2 ) + 4s 5 (t 2 + tu + u 2 ) + s 6 (t + u) + 

+ 2stw(t 2 + iu + m 2 ) 2 + s 2 (t + w) (t 4 + At 3 u + 3t 2 u 2 + Atu 3 + m 4 ) 

+ 2s 3 (2t 4 + 5£ 3 m + 5tV + 5tu 3 + 2m 4 ) + 2s 4 (t + u) (?,t 2 + 2tu + 3m 2 ) 

+ 4s 5 (t 2 + tu + m 2 ) + s 6 (t + M 



t 2 u 2 (t + uY 



da / 



lo(8) 







3e^ 4 O 



lo(8) 



Q2« 



327rMst(s + t) 2 (s + m) 2 (£ + m) 2 



x 4u 4 + 8tM 3 + 13tV + 10t 3 M + 5t 4 + 2s(t + uf (5t + 4m) + 13s 2 (t + uf 



+ 8s 3 (* + «) + 4s 4 + 
+ 8s 3 (t + u) + 4s 4 



4tM (2t 2 + 3tu + 2m 2 ) + 8s(t + m) 3 + 12s 2 (t + uf 



da / 



9t (O 



3 S 



,(8) 



g 2 M 



36ns 2 (s + t) 2 (s + u) 2 (t + u) 2 



X 



tV + StM (t + u) + S 2 (t 2 + tu + M 2 ) + £ a £ t 



— t M — StM (£ + m) + S t« 



da 
dt 



5e 2 gUo M 8) ])Q 



(A.4) 



(A.5) 



(A.6) 



a u A + 4tM 5 + M f 



( 7 ( ? ^ CC [P 1 ](?) 97rMs2(s + t)3(s + M)3(t + M) 3 
X \^tu{t + M) 3 (t 2 + tu + M 2 ) + S (t 6 + 4t 5 M + lltV + 14*V + lit 

+ s 2 (t + m) (4t 4 + 7t 3 M + 9t 2 u 2 + 7tM 3 + 4m 4 ) 

+ s 3 (7t 4 + 14t 3 M + 16tV + 14tM 3 + 7m 4 ) + s 4 {t + u) (7t 2 + Atu + 7m 2 ) 

+ 4s 5 (t 2 + tu + m 2 ) + s 6 (t + u) + ^ b t 2 u 2 {t + m) 3 + 2stM(t 2 + tu + m 2 ) 2 

+ s 2 (t + u) (t 4 + 4t 3 M + 3t V + 4tM 3 + m 4 ) + 2s 3 (2t 4 + 5t 3 M + 5fV + 5£m 3 + 2m 4 ) 

+ 2s 4 (t + u) (3t 2 + 2tu + 3m 2 ) + 4s 5 (t 2 + tu + m 2 ) + s 6 (t + 



(A.7) 



da 
dt 

x 



3p(8) 



( [3 (8)1 \ 

(19^CC[P J^J- 87rM 3 st(s + t) 4 (s + M) 4 (t + M) 4 

t 2 (t + m) 2 (4t 6 + 20t 5 M + 105t 4 M 2 + 142t 3 M 3 + 125* V + 72tM 5 + 36m 6 ) 
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+ 2st (t + u) (lAt 7 + 143t 6 M + 463fV + 707tV + 679t 3 M 4 + 426tV + 180tM 6 + 36m 7 ) 
+ s 2 (l49t 8 + 1212t 7 u + 3692tV + 6228tV + 6597tV + 4424£ 3 m 5 + 1838t 2 M 6 



+ 432tw 7 + 36m 8 



+ 4s 3 (t + u) (93t 6 + 492t 5 M + 1065t 4 M 2 + 1293t 3 M 3 + 818t 2 u 4 + 267£m 5 + 36m 6 ) 
+ s 4 (514t 6 + 2772t 5 M + 6597* V + 8444tV + 5816t 2 M 4 + 2040tM 5 + 297m 6 ) 
+ 2s 5 (t + u) (232t 4 + 873t 3 M + 1339tV + 83Km 3 + 189m 4 ) 

+ s 6 (305t 4 + 1212t 3 M + 1838tV + 1212tM 3 + 297m 4 ) + 144s 7 (t + m) 3 + 36s 8 (t + u) 

+ - 4t 3 M 2 (t + m) 2 (8t 3 - 8t 2 M - 25tw 2 - 18m 3 ) 

- 8st 2 M (t + u) (8t 5 + 4t 4 M - 72tV - 140tV - 105tM 4 - 27m 5 ) 
+ 8s 2 t (-At 7 - 12t 6 M + 103tV + 440tV + 687tV + 535tV + 202tM 6 + 27m 7 
+ 8s 3 (-4t 7 + 68t 6 M + 434tV + 923tV + 970t 3 M 4 + 532tV + 132£m 6 + 9m 7 ) 
+ 4s 4 (33t 6 + 412t 5 M + 1303iV + 1858t 3 M 3 + 1377tV + 498£m 5 + 63m 6 ) 
+ 8s 5 (t + u) (38t 4 + 198t 3 M + 302tV + 201tM 3 + 45m 4 ) 

+ 8s 6 (35t 4 + 147t 3 M + 221tV + 147£m 3 + 36m 4 ) + 144s 7 (t + m) 3 + 36s 8 (t + u) 



(A.8) 



da 
~dt 



(l9 



cc 



3p(8) 



e 2 gUoM 8) ])Q 



8nM 3 s 2 (s + t)\s + u)\t + u) A 



x \t 4 u 2 {t + m) 3 (t 2 + 2tu + 5m 2 ) + st 3 u 2 (t + m) 2 (3t 3 + 18t 2 M + 47£m 2 + 12m 3 

+ s 2 t 2 (t + u) (t 6 + 2t 5 M + lOtV + 82t 3 M 3 + 183tV + lOOtM 5 + 22m 6 ) 

+ sh (5t 7 + 24t 6 M + 92tV + 296tV + 465tV + 316tV + 94£m 6 + 12m 7 ) 

+ s 4 (t + m) (l4t 6 + 72t 5 M + 193t 4 M 2 + 272tV + 116tV + 16tM 5 + 5m 6 ) 

+ s 5 (22t 6 + 124t 5 M + 283tV + 316tV + 132tV + 8tM 5 - 5m 6 ) 

+ s 6 (t + m) 2 (l7t 3 + 37t 2 M + 31tM 2 - 5m 3 ) + s 7 (5t 4 + 12t 3 u + 22t 2 M 2 + 12£m 3 + 5m 4 ) 

+ Z a Z b - 4t 4 u 4 (t + m) 3 - 4stV(t + m) 2 (2t 2 + 9tu + 2m 2 ) 

- 4s 2 tV (t + u) (t 4 + 14t 3 M + 34tV + 16tM 3 + 3m 4 ) 

- 4s 3 tM 2 (3t 5 + 28t 4 M + 56tV + 38tV + 11£m 4 + 2m 5 ) 
+ 4s 4 tM (t + m) (2t 4 + 5t 3 M - 8t V - 5tM 3 - 2m 4 ) 
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+ As 5 t 2 u (Qt 3 + 21t 2 u + 20tu 2 + 9m 3 ) + 8s 6 tu(t + uf (3t + u) 
+ 8s 7 tu (t 2 + tu + m 2 ) 



(A.9) 



3p(8) 



= 



-e 2 <? 4 (o 



SpW 



<5c 



dcr / 

x | - 3t 5 M 2 (t + m) 3 (t 2 + 2tu + 5m 2 ) 

- st 2 u(t + u) 2 (8t 6 + 49t 5 u + 186tV + 317tV + 256tV + lUtu 5 + 48m 6 ) 

- s 2 t{t + u) 2 (3t 7 + 59t 6 M + 387tV + 979tV + 1090tV + 738tV + 384tw 6 + 96m 7 ) 
+ s 3 (-15t 9 - 292t 8 M - 1812tV - 4840t 6 M 3 - 701KV - 6384t 4 M 5 - 4090t 3 M 6 

- 1900tV - 528tM 8 - 48m 9 ) 

- s 4 (t + m) (42t 7 + 696t 6 M + 2739t 5 M 2 + 4272tV + 3804tV + 2368tV + 1071£m 6 
+ 192m 7 ) 

+ s 5 (-66t 7 - 1076t 6 M - 3897*V - 6384tV - 6172tV - 4144tV - 1881£m 6 

- 396m 7 ) 

- s 6 (t + u) (51t 5 + 922t 4 M + 2028tV + 2062tV + 1377£m 4 + 504m 5 ) 
+ s 7 (-15t 5 - 592t 4 M - 1602t 3 M 2 - 1900tV - 1263tM 4 - 396m 5 ) 

- 48s 8 m (t + u) (5t 2 + 7tM + 4m 2 ) - 48s 9 M(t + m) 2 + £ a £ h 

+ 4st 3 M 3 (t + m) 2 (2t + 3m) (5t 2 - 2tu - 4m 2 ) 
+ 4s 2 t 3 M 2 (t + m) 2 (lit 3 + 47t 2 M + 3Hm 2 - 3m 3 ) 



12*V(* + u) i 



4s 3 tM (At' + 81t D M + 284t°M z + 472t 4 M d + 486t<V + 341t z u° + 158£m d + 36m 



+ 4s 4 m (t + u) (34t 6 + 225t 5 M + 424tV + 495fV + 386tV + 180tM 5 + 24m 6 ) 
+ 4s 5 m (l02t 6 + 469t 5 M + 892tV + 1033t 3 M 3 + 828t 2 M 4 + 408£m 5 + 84m 6 ) 
+ 8s 6 m (t + u) (77t 4 + 186t 3 M + 209tV + 156tM 3 + 60m 4 ) 
+ 8s 7 m (65t 4 + 183t 3 M + 221tV + 150£m 3 + 48m 4 ) 

+ 48s 8 m (t + u) (5t 2 + 7tu + 4m 2 ) + 48s 9 m(£ + uf 



(A.10) 



dcr / 



lo(8) 



Q 



12hMsH(s + uY 
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X lu 2 + S 2 + ^ 



-u 2 + s 2 



(A.11) 



cc 



3 ^ 8) ] ?) 



° 2 5s 4 (c[ 3 ^ 8) ])g 



x {u 2 + 2tu + It 2 + 2st + s 2 + £ a £ b 



da / 
- ( 7 g - cc 



3d(8) 



72itM 3 s 3 u 

u(2t + u) -2st-s 2 j. 

-e 2 W(9 [[ 3 P (8) l)Q 2 (3s + 2t + 3M) s 



(A.12) 



x ^ m 2 + s 2 + 



9 = 



2 i 2 
— M + S 



9ttM 3 s 2 £(s + uy 



(A.13) 



d(7 

~dt 



(yq -> cc ppf 85 ] g) 



9ttM 3 s 2 (s + m) 4 



x tu 2 + 4sm (t + u) + s 2 (t + 4m) + £ a £ 6 



tM 2 + S 2 t 



(A.14) 



da 



3p(8) 
-<2 



= 



2e 2 W0 



g 



x j - m 2 (7t 2 + 12*it + 6m 2 ) - 12su(f + m) 2 + s 2 (-7t 2 - 24£m - 12m 2 ) - 12s 3 (t + u) 
- u 2 (ht 2 + 12tu + 6m 2 ) - 12sm 2 (t + u) + s 2 t (5t + 12m) + 12s 3 (t + u) 

(A.15) 



+ 6s 
da 



(99 



cc 



1 Q (1) 



9U0 



<~>n 



647rMs 3 t(s + t) 2 M(s + m) 2 (* + m) 2 



x jtV (4t 4 + 8t 3 M + 13tV + 8tM 3 + 4m 4 ) 

+ 2stw (t + u) (4t 4 + 12t 3 M + 21iV + 12tM 3 + Au A 
+ s 2 (4t 6 + 32t 5 M + llltV + 168t 3 M 3 + 11 ltV + 32tM 5 + 4m 6 ) 
+ 2s 3 (t + u) (At 4 + 29t 3 M + 55tV + 29£m 3 + 4m 4 ) 

+ s 4 (l3t 4 + 66t 3 M + llltV + 66tM 3 + 13m 4 ) + 8s 5 (t + u) (t 2 + 3tM + m 2 ) 
+ 4s 6 (t + m) 2 + i a i b 4tV (2t 2 + 3tu + 2m 2 ) + 8stV (t + u) (3t 2 + 5tu + 3m 2 ) 
+ 8s 2 tM(t + m) 2 (3t 2 + 7tu + 3m 2 ) + 8s 3 (t + u) (t 4 + 7t 3 u + 13tV + 7£m 3 + m 4 ) 

+ 4s 4 (t 2 + 3tu + m 2 ) (3t 2 + 7tM + 3m 2 ) + 8s 5 (f + u) (t 2 + 3tM + m 2 ) + 4s 6 (t + uf 
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da / 



cc 



M 



2lQirs 2 (s + t) 2 (s + u) 2 (t + u) 2 



x 



da 
dt 



t 2 U 2 + StM (* + «) + S 2 {f +tu + M 2 ) + £ a £ fe 



t u — stu (t + u) + s tu 



>, (A.17) 



X 



j*u(i + m) 3 (t 2 + + u 2 ) + s (t 6 + 4t 5 u + lltV + Ut 3 u 3 + lltV + 4tu 5 + m 6 ) 
+ s 2 (t + u) (4t 4 + 7t 3 u + 9tV + 7tw 3 + 4m 4 ) 

+ s 3 (7t 4 + 14£ 3 m + 16tV + Utu 3 + 7m 4 ) + s 4 (t + u) (7t 2 + 4tM + 7m 2 ) 
+ 4s 5 (t 2 + tu + m 2 ) + s 6 (t + u) + £ a £ 6 t 2 u 2 {t + u) 3 + 2stu(t 2 + tu + m 2 ) 2 

+ S 2 (t + u) (t 4 + 4t 3 M + 3tV + 4tM 3 + M 4 ) 

+ 2s 3 (2t 4 + 5t 3 u + 5t V + 5tu 3 + 2m 4 ) + 2s 4 (t + u) (3t 2 + 2tu + 3m 2 ) 
+ 4s 5 (t 2 + tw + w 2 ) +s 6 (t + w) |, 

9 6 s (0 



cc 



3p(l) 



123tV 



56£m 5 + 12w 6 ) 



da / 

— [gg .^ ro Jyy - 487rM 3 s 3 t ( s + t )4 M(s + M)4(t + M) 4 

x |9t 4 M 4 (t + m) 2 (4t 4 + 8t 3 M + 13tV + 8tM 3 + 4m 4 ) 

+ 12st 3 M 3 (t + m) (l2t 6 + 56t 5 M + 123tV + 163t 3 M 3 
+ 4s 2 tV (55t 8 + 439t 7 M + 1517t 6 M 2 + 3076t 5 M 3 + 3891fV + 3076tV 
+ 1517t 2 M 6 + 439tw 7 + 55m 8 ) 

+ 2s 3 tw (t + u) (j2t 8 + 806t 7 M + 3419tV + 7828fV + 10398tV + 7828tV 
+ 3419t 2 M 6 + 806tw 7 + 72m 8 ) 

+ s 4 (36t 10 + 816t\ + 6068tV + 22494tV + 48086tV + 61808t 5 M 5 
+ 48086tV + 22494tV + 6068tV + 81Qtu 9 + 36m 10 ) 
+ 4s 5 (t + u) (36t 8 + 501t 7 M + 2575fV + 6538t 5 M 3 + 8914tV + 6538fV 
+ 2575tV + 50Hm 7 + 36m 8 ) 

+ s 6 (297t 8 + 3432t 7 M + 15564fV + 36452tV + 48086tV + 36452tV 
+ 15564t 2 M 6 + 3432tM 7 + 297m 8 ) 



(A.18) 
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+ 2s 7 (t + u) (l89t 6 + 1527t 5 M + 4625t 4 M 2 + 6622tV + 4625£ 2 m 4 + 1527£m 5 + 189m 6 ) 
+ s 8 (297t 6 + 2U8t 5 u + 6068t 4 M 2 + 8450tV + 6068tV + 2148£m 5 + 297m 6 ) 



+ 4s y (t + u) (36t 4 + \m 6 u + 271fV + 168£m 3 + 36m 
+ 4s 10 (9t 4 + 36£ 3 m + 55tV + 36tu 3 + 9m 4 ) + 36tV(t + m) 2 (2t 2 + 3tu + 2m 2 ) 
+ 48st 4 M 4 (t + m) (7t 4 + 27t 3 M + 39tV + 27£m 3 + 7m 4 ) 

+ 4sVm 3 (l62t 6 + 1057t 5 M + 2712tV + 3630t 3 M 3 + 2712tV + 1057tM 5 + 162m 6 ) 
+ 8s 3 tV (t + m) (81t 6 + 631t 5 M + 1772tV + 2424iV + 1772tV + 631£m 5 + 81m 6 ) 
+ 4s 4 tw (84t 8 + 1057t 7 w + 4806t 6 M 2 + 10979t 5 M 3 + 14284* V + 1097&V 
+ 4806t 2 w 6 + 1057tw 7 + 84m 8 ) 



+ 8s 5 {t + u) [9t 8 + 195t Y M + 1161* V + 3041tV + 4118* V + 3041*V 



+ 1161* 2 m 6 + 195*m 7 + 9m 8 
+ 4s 6 (63t 8 + 792* 7 m + 3672* V + 8570* 5 m 3 + 11250*V + 8570* V 
+ 3672*V + 792*m 7 + 63m 8 ) 

+ 24s 7 (* + m) 3 (l5t 4 + 93* 3 m + 169*V + 93*m 3 + 15m 4 ) 
+ 4s 8 (* + uf (72t 4 + 384* 3 m + 641*V + 384*m 3 + 72m 4 ) 

+ 48s 9 (* + m) 3 (3* 2 + 8*m + 3m 2 ) + 36s 10 (* + uf 
da, _ r ,„ mn x 9 6 s (O 



(A.19) 



cc 



spW 



48ttM 3 s 2 (s + tf(s + u)\t + m) 4 



X j*V (* + U) {f + 4* 5 M + 12*V - 2*V + 12*V + 4tM 5 + M 6 ) 

+ s*V (3t 6 + 32* 5 m + 58* 4 m 2 + 48* 3 m 3 + 58* V + 32*m 5 + 3m 6 ) 
+ s 2 (* + m) (* 8 + 2* 7 m + 22*V + 46*V + 10*V + 46*V + 22*V + 2*m 7 + m 8 ) 
+ s 3 (5* 8 + 32* 7 m + 68*V + 16*V - 74*V + 16*V + 68* V + 32*m 7 + 5m 8 ) 
+ 2s 4 (* + u) (8* 6 + 21* 5 m + 7* 4 m 2 - 44* 3 m 3 + 7* 2 m 4 + 21*m 5 + 8m 6 ) 



2s s ( 5* b + 24*°m + 28* 4 m 2 + 8* a M 3 + 28* z m 4 + 24*m° + 5m 



3„,3 



+ 2s 6 (* + u) (t 2 + m 2 ) (5* 2 + 24*m + 5m 2 ) + 8s 7 (2* 4 + 4* 3 m + 3* V + 4*m 3 + 2m 4 ) 
+ s 8 (* + u) (5* 2 - 2tu + 5m 2 ) + s 9 (* 2 + m 2 ) + - 4* 4 (* - m) 2 m 4 (* + u) 
- Ast 3 u 3 (2t 4 + t 3 M - 4t 2 M 2 + tM 3 + 2m 4 ) - 4s 2 tV (t 5 + 3t 4 M + 3tM 4 + m 5 ) 
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+ 4s 3 tV(t 2 + m 2 ) 2 + AsHu (t + u) (2t A + t 3 u + 2t 2 u 2 + tu 3 + 2m 4 ) 
- 4s 5 t V {f + u 2 ) - BsHu (t + u) (f + u 2 ) 



(A.20) 



do" / _ r 3 „ri)i \ (C 



3p(l) 



2407rM 3 sH(s + t) 4 w(s + uf(t + m) 4 



x | - 12tV(t + uf (4t 4 + 8t 3 M + 13fV + 8£m 3 + 4m 4 ) 

- 3stV (t + u) (65t 6 + 292t 5 M + 572tV + 710t 3 M 3 + 572t 2 M 4 + 292£m 5 + 65m 6 ) 

- sH 2 u 2 (296t 8 + 2273t 7 M + 6832tV + 11846tV + 13992£ 4 m 4 + 11846t 3 M 5 
+ 6832t 2 M 6 + 2273£m 7 + 296m 8 ) 

- sHu (t + u) (l95t 8 + 2078t 7 M + 7202£ 6 m 2 + 12562tV + 14718tV 
+ 12562t 3 M 5 + 7202tV + 2078tM 7 + 195m 8 ) 

+ s A (-48t 10 - 107K 9 m - 6832tV - 19764tV - 33496tV - 39106tV 

- 33496£ 4 m 6 - 19764iV - 6832tV - 107Km 9 - 48m 10 ) 

- 2s 5 (t + u) (96t 8 + 1200t 7 M + 4723tV + 8917t 5 M 3 + 10636tV + 8917t 3 M 5 
+ 4723tV + 1200tM 7 + 96m 8 ) 

- 2s 6 (l98t 8 + 1923t 7 M + 6996tV + 13640tV + 16748* V + 13640tV 



+ 6996t 2 M 6 + 1923tM 7 + 198m 8 



- 2s 7 (t + m) (3t 2 + Atu + 3m 2 ) (84t 4 + 445t 3 w + 740* V + 445*m 3 + 84m 4 ) 

- 4s 8 (99t 6 + 648t 5 M + 1708t 4 M 2 + 2320t 3 M 3 + 1708* 2 m 4 + 648*m 5 + 99m 6 ) 

- s 9 (t + u) (l92t 4 + 879t 3 M + 1394*V + 879*m 3 + 192m 4 ) 
+ s 10 (-48t 4 - 195t 3 M - 296*V - 195*m 3 - 48m 4 ) 



+ £a£b -48tV(t + M) 2 (2t 2 + 3tM + 2M 2 ) 

- 12stV (t + u) (24t 4 + 103t 3 M + 154*V + 103*m 3 + 24m 4 ) 

- 4s 2 tV (48t 6 + 506t 5 M + 1557fV + 2196t 3 M 3 + 1557*V + 506*m 5 + 48 
+ 4s 3 tV (t + u) (48t 6 - 37t 5 M - 794t 4 M 2 - 1374tV - 794*V - 37*m 5 + 48m 
+ 4s 4 tM (72t 8 + 520t 7 M + 1077t 6 M 2 + 716tV + 178fV + 716t 3 M 5 + 1077* 
+ 520*m 7 + 72m 8 ) 

+ 4s 5 (* + u) (24t 8 + 360* 7 m + 1314*V + 1957*V + 1954*V + 1957*V 



^m 6 



2 m 6 



18 



+ 1314tV + 360£m 7 + 24m 8 ) 

+ 4s 6 (84t 8 + 792t 7 u + 2688t 6 M 2 + 4811£ 5 m 3 + 5664tV + 4811t 3 M 5 + 2688tV 



792tu 7 + 84m 8 



+ 24s 7 (t + u) (20t 6 + 128t 5 M + 309t 4 M 2 + 396t 3 M 3 + 309tV + 128tM 5 + 20m 6 ) 
+ 64s 8 (t + uf (6t 4 + 27t 3 M + 41tV + 27£m 3 + 6m 4 ) 

+ 96s 9 (t + uf (2t + u)(t + 2u) + 48s 10 (t + u)^ 



(A.21) 



da 
~dt 



(gg^cc[S Q ]g) \02A>kMsH(s + t) 2 u(s + u) 2 (t + u) 2 

x jiV (4t 4 + 8t 3 u + 13tV + 8tu 3 + 4m 4 ) 

+ 2stu (t + u) (2t 4 + 8t 3 w + 15* V + 8tw 3 + 2w 4 ) 

+ s 2 (4t 6 + 20t 5 M + 71tV + 108t 3 M 3 + 71tV + 20tM 5 + 4m 6 ) 



+ 2s 3 (t + u) (4t 4 + 19t 3 M + 35iV + 19tM 3 + 4m 
+ s 4 (l3t 4 + 46t 3 M + 71tV + 46tw 3 + 13m 4 ) + 4s 5 (t + u) (2t 2 + 3tM + 2m 2 ) 

+ 4s 6 (t 2 + tu + m 2 ) + £ a £ b At 3 u 3 (2t 2 + 3tu + 2m 2 ) + 4si V (t + m) (4t 2 + 7tu + 4m 2 ) 

+ 4s 2 tM (4t 4 + 16t 3 M + 25t V + 16tM 3 + 4m 4 ) 
+ 4s 3 (t + u) (2t 4 + 9t 3 M + 16t V + 9tM 3 + 2m 4 ) 

+ 4s 4 (3t 4 + 1H 3 m + 17tV + lltM 3 + 3m 4 ) + 4s 5 (t + u) (2t 2 + 3tu + 2m 2 ) 

(A.22) 

35(8)" 



u 



+ 4s 6 (t 2 + tu 

3 *i 8) ] 9) = 



da / 



cc 



1474567rM 5 s 3 t(s + t) 2 M(s + u) 2 (t + m) 2 



x |27t 2 M 2 (t + uf (6t 4 + 9t 3 M + 2tV + 9tM 3 + 6m 4 ) 

- stu (t + u) (l08t 6 - 297t 5 M - 4349tV - 3352tV - 4349tV - 297£m 5 + 108m 6 ) 
+ s 2 (l62t 8 + 189t 7 M + 3728tV + 13757t 5 M 3 + 15792tV + 13757tV + 3728fV 
+ 189tM 7 + 162m 8 ) 

+ s 3 (t + m) (567t 6 + 4079t 5 M + 9678iV + 11636tV + 9678* V + 4079*m 5 + 567m 6 ) 
+ s 4 (702t 6 + 7701t 5 M + 15792t 4 M 2 + 21314* 3 m 3 + 15792tV + 7701*m 5 + 702m 6 ) 
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+ s 5 (t + u) (594t 4 + 7107t 3 u + 6650tV + 7107£m 3 + 594m 4 ) 

+ s 6 (702t 4 + 4646t 3 M + 3728tV + 4646£m 3 + 702m 4 ) + 189s 7 (t + w) (3t 2 - 2tM + 3m 2 ) 
+ 54s 8 (3t 2 - 2£m + 3w 2 ) + £ a £ b - 128st 3 M 3 (t + u) (l9t 2 + Utu + 19m 2 ) 

- 128s 2 tV (l9t 4 + 79t 3 M + 93tV + 79tM 3 + 19m 4 ) 

- 128s 3 tV (t + m) (lit 2 + 52tM + 11m 2 ) 



- 5888s 4 t 3 M 3 + 1408s & tV (t + u) + 2432s 6 £ V 



(A.23) 



A (8) ] g) 



gUO ^ 



1 D (8) 



da / 

— [gg -> cc L , 1 j y; - 23 047rM 3 s 3 t ( s + t )3^ s - M ^ 3 ^ - 
x |27t'V (i + m) (4t 4 + 8t 3 M + 13tV + 8tM 3 + 4m 4 ) 

+ st 2 u 2 (t + m) 2 (l52t 4 + 780t 3 M + 1229tV + 780tM 3 + 152m 4 ) 
+ s 2 tM (t + u) (l52t 6 + 328t 5 M + 2479tV + 3735t 3 M 3 + 2479tV + 328tM 5 + 152m 6 ) 
+ s 3 (l08t 8 + 1084t 7 M + 2807tV + 6262t 5 M 3 + 7694tV + 6262t 3 M 5 + 2807tV 
+ 1084tM 7 + 108m 8 ) 

+ s 4 {t + u) (324t 6 + 2617t 5 M + 3597tV + 4097*V + 3597tV + 2617tM 5 + 324m 6 ) 
+ s 5 (567t 6 + 4018t 5 M + 6214tV + 6262tV + 6214* V + 4018£m 5 + 567m 6 ) 
+ s 6 (t + u) (567t 4 + 2374t 3 M + 433fV + 2374£m 3 + 567m 4 ) 

+ 4s 7 (81t 4 + 271t 3 M + 120tV + 271tM 3 + 81m 4 ) + 4s 8 (t + u) (27t 2 + Utu + 27m 2 ) 

+ CaCb 108t 4 M 4 (t + u) (2t 2 + 3tu + 2m 2 ) + 4st 3 M 3 (t + uf (ll9t 2 + 265£m + 119m 2 ) 

+ 8s 2 tV (t + u) (65t 4 + 176t 3 M + 34H 2 m 2 + 176£m 3 + 65m 4 ) 

+ 4s 3 tM (ll9t 6 + 520t 5 M + 824tV + 1138£ 3 m 3 + 824* V + 520£m 5 + 119m 6 ) 

+ 4s 4 (t + u) (54t 6 + 487t 5 M + 667tV + 602t 3 M 3 + 667tV + 487tM 5 + 54m 6 ) 

+ 4s 5 (l35t 6 + 866t 5 M + 1312tV + 1194tV + 1312tV + 866£m 5 + 135m 6 ) 

+ 4s 6 (t + u) (l35t 4 + 542t 3 M + 88iV + 542£m 3 + 135m 4 ) 

+ 4s 7 (81t 4 + 271t 3 M + 120tV + 27Hm 3 + 81m 4 ) 

+ 4s 8 (t + M) (27t 2 + 1Um + 27m 2 ) j, (A.24) 

da, _ r ,„ mn \ S^Ofpf 

dt (" 



cc 



3p(8) 



768vrM 3 s 3 t(s + tfu(s + uf(t + uf 
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x jgtV(t + uf (At 4 + 8t 3 u + 13* V + 8tu 3 + 4m 4 ) 

+ 12st 3 u 3 (t + u) (9t 6 + 43t 5 M + 97t 4 u 2 + 131tV + 97*V + 43*m 5 + 9m 6 ) 
+ 4s 2 tV (37t 8 + 295* 7 m + 1045* 6 m 2 + 2185f 5 u 3 + 2801t 4 M 4 + 2185tV 
+ 1045tV + 295*m 7 + 37m 8 ) 

+ 2s 3 tu (t + u) (54t 8 + 536t 7 M + 2261tV + 5324*V + 7178t 4 u 4 + 5324* 
+ 2261* 2 m 6 + 536*m 7 + 54m 8 ) 

+ 2s 4 (l8* 10 + 312* 9 m + 2090*V + 7585*V + 16243*V + 20924* 5 m 
+ 16243* 4 m 6 + 7585*V + 2090* 2 m 8 + 312*m 9 + 18m 10 ) 
+ 4s 5 (* + u) (36t 8 + 384* 7 m + 1801*V + 4450* 5 m 3 + 6012* V + 4450*V 
+ 1801*V + 384*m 7 + 36m 8 ) 

+ s 6 (297t 8 + 2736* 7 m + 11204* 6 m 2 + 25004*V + 32486*V + 25004*V 
+ 11204* 2 m 6 + 2736*m 7 + 297m 8 ) 

+ 2s 7 (* + m) (l89t 6 + 1179* 5 m + 3191*V + 4394*V + 3191*V + 1179*m 5 + 189m 
+ s 8 (t + m) 2 (297t 4 + 1086* 3 m + 1711*V + 1086*m 3 + 297m 4 ) 
+ 4s 9 (* + u) (36t 4 + 120* 3 m + 175*V + 120*m 3 + 36m 4 ) 

+ 4s 10 (9t 4 + 27t 3 M + 37*V + 27*m 3 + 9m 4 ) + ^6 36*V(* + m) 2 (2* 2 + 3*m + 2m 2 ) 

+ 12s* 4 m 4 (f + u) (23t 4 + 89t 3 M + 128*V + 89*m 3 + 23m 4 ) 
+ 4s 2 t 3 M 3 (ll7t 6 + 755* 5 m + 1920* 4 m 2 + 2560*V + 1920*V + 755*m 5 + 117m 6 ) 
+ 4s 3 *V (* + m) (ll7t 6 + 851t 5 M + 2312*V + 3128fV + 2312*V + 851*m 5 + 117m 6 ) 
+ 4s 4 *m (69t 8 + 759* 7 m + 3185*V + 6993*V + 8992*V + 6993* 3 m 
+ 3185* 2 m 6 + 759*m 7 + 69m 8 ) 

+ 4s 5 (* + u) (l8t 8 + 318* 7 m + 1618*V + 3911* 5 m 3 + 5198*V + 3911* 
+ 1618* 2 m 6 + 318*m 7 + 18m 8 ) 

+ 4s 6 (63t 8 + 645* 7 m + 2592*V + 5631* 5 m 3 + 7246* 4 m 4 + 5631fV 
+ 2592* V + 645*m 7 + 63m 8 ) 

+ 4s 7 (* + m) (90t 6 + 570* 5 m + 1481*V + 2030*V + 1481*V + 570*m 5 + 90m 6 
4s 8 (72t 6 + 411* 5 m + 1005* 4 m 2 + 1336*V + 1005* 2 m 4 + 411*m 5 + 72m 6 ) 



*3..5 



+ 
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+ 48s 9 (t + m) 3 (3t 2 + Atu + 3m 2 ) + 36s 10 (t + uf (t 2 + tu + m 2 ) 

5g* (O [«P< 8 >]) 



(A.25) 



da / 
dt 



3p(8) 



768ttM 3 s 2 (s + t) 4 (s + u)\t + m) 4 



x jtV (t + u) (t 6 + 4t 5 M + 16tV + 6t 3 M 3 + 16tV + Atu b + u 
+ st 2 u 2 (t 2 + tu + m 2 ) (3t 4 + 37t 3 M + 78tV + 37tM 3 + 3m 4 ) 

+ s 2 (t + u) (t 8 + 2t 7 u + 38tV + 150tV + 218tV + 150tV + 38tV + 2tu 7 + m 8 ) 

+ s 3 (5t 8 + 40t 7 M + 188t 6 M 2 + 472tV + 622tV + 472t 3 M 5 + 188tV + AOtu 7 + 5m 8 ) 

+ 2s 4 (t + u) (lOt 6 + 49t 5 M + 135tV + 176t 3 M 3 + 135tV + 49tM 5 + 10m 6 ) 

+ 2s 5 (lit 6 + 76t 5 w + 184tV + 236tV + 184tV + 76tw 5 + 11m 6 ) 

+ 2s 6 (t + u) (lit 4 + 48t 3 M + 46tV + 48tw 3 + 11m 4 ) 

+ 20s 7 (t + u) 2 (t 2 + m 2 ) + s 8 (t + u) (5t 2 - 2tM + 5m 2 ) + s 9 (t 2 + m 2 ) 

+ iaib - 8tV (t + u) (t 2 + m 2 ) - 8st 3 w 3 (2t 4 + 8t 3 M + lltV + 8tM 3 + 2m 4 ) 

- 16s 2 tV (t + u) (t 2 + 3tu + w 2 ) 2 

- 8s 3 tw (t 6 + 8t 5 M + 31tV + 47tV + 31tV + 8tw 5 + m 6 ) 

- 8s 4 tM (t + u) (t 4 + 6t 3 M + 12tV + 6tw 3 + m 4 ) + 8s 5 tV (3t 2 + 5tM + 3m 2 ) 

+ 8s 6 tM (t + u) (t 2 + 4tM + m 2 ) + 8s 7 tM (t 2 + tu + m 2 ) \ , (A.26) 



da / 



cc 



3p(8) 
*1 



-9UO 



Sp(8) 



768vrM 3 s 3 t(s + t) 4 M(s + M) 4 (t + m) 4 



x I - 12t 4 M 4 (t + m) 2 (4t 4 + 8t 3 M + 13t 2 M 2 + 8tM 3 + 4m 4 ) 

- 3st 3 M 3 (t + m) (49t 6 + 220t 5 M + 432t 4 M 2 + 542t 3 M 3 + 432t 2 M 4 + 220tM 5 + 49m 6 ) 

- s 2 tV (200t 8 + 1577t 7 M + 4808t 6 M 2 + 8498t 5 M 3 + 10144tV + 8498t 3 M 5 
+ 4808t 2 M 6 + 1577tM 7 + 200m 8 ) 

- s 3 tM (t + m) (l47t 8 + 1430t 7 M + 5066t 6 M 2 + 9242tV + 11150t 4 M 4 



+ 9242t 3 M 5 + 5066t 2 M 6 + 1430tM 7 + 147m 8 



+ s 4 (-48t 10 - 807t 9 M - 4808t 8 M 2 - 14308t 7 M 3 - 25552tV - 30602t 5 M 5 
- 25552t 4 M 6 - 14308t 3 M 7 - 4808t 2 M 8 - 807tM 9 - 48m 10 ) 
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- 2s 5 (f + u) (96t 8 + 882t 7 M + 3367tV + 6829tV + 8472tV + 6829tV 
+ 3367t 2 M 6 + 882tw 7 + 96m 8 ) 

- 2s 6 (l98t 8 + 1461t 7 M + 5072t 6 M 2 + 10196tV + 12776* V + 10196£ 3 m 5 
+ 5072tV + 146Km 7 + 198m 8 ) 

- 2s 7 (t + u) (252t 6 + 1209t 5 M + 3040* V + 4114t 3 M 3 + 3040tV + 1209£m 5 + 252m 6 ) 

- 4s 8 (t + uf (99t 4 + 291t 3 M + 521tV + 291£m 3 + 99m 4 ) 



- s 9 (t + u) (I92t 4 + 615t 3 M + 962tV + 615£m 3 + 192m 
+ s 10 (-48t 4 - 147t 3 M - 200tV - 147tw 3 - 48m 4 ) 



- 48tV(t + uf (2t 2 + 3tu + 2m 2 ) 



- 24stV (t + m) (3t 2 + 5tu + 3m 2 ) (3t 



2 + 8tu + 3u 2 



8s 2 tV ( 15t b + 170t°w + 525tV 4 + 739fV + 525t z M 4 + 170tM & + 15m 



*4„.2 



3„.3 



+ 8s 3 tV (t + u) (l5t 6 - 5t 5 M - 242tV - 422tV - 242tV - 5£m 5 + 15m 6 ) 
+ 8s 4 tM (3t 2 + 7tu + 3m 2 ) (9t 6 + 40t 5 M + 30tV - t 3 u 3 + 30tV + 40£m 5 + 9m 6 ) 
+ 8s 5 (t + u) (l2t 8 + 135t 7 M + 455tV + 688t 5 M 3 + 694fV + 688t 3 M 5 
+ 455tV + 135tM 7 + 12m 8 ) 

+ 8s 6 (42t 8 + 300t 7 M + 924iV + 1617t 5 M 3 + 1901tV + 1617fV 
+ 924t 2 M 6 + 300tM 7 + 42m 8 ) 

+ 16s 7 (t + u) (30t 6 + 135t 5 M + 308tV + 395tV + 308tV + 135£m 5 + 30m 6 ) 
+ 8s 8 (48t 6 + 231t 5 M + 549tV + 731fV + 549tV + 231£m 5 + 48m 6 ) 

+ 24s 9 (t + u) (8t 4 + 25t 3 M + 38tV + 25tM 3 + 8m 4 ) + 48s 10 (t + m) 2 (t 2 + tu + m 2 ) 



da 
~dt 



(gq -> cc 



<~>n 



72%Ms 2 t(s + uf 



X{u 2 + S 2 + CaCb 



2 , 2 
U + S 



(A.27) 



(A.28) 



da / 



3 D (1) 



-9t O 



3p(l)' 
-in 



(3s + 2t + 3uf 



547rM 3 s 2 t(s + m) 



x {u 2 + s 2 + CaCl 



2 i 2 
— M + S 



(A.29) 
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da 
dt 



(gq -> cc [ 3 P! (1) ] q) 



-9! 



o[-p«]) 



27ttM 3 s 2 (s + m) 



x < tM 2 4" 4sm (t + u) + s 2 (t + 4m) + £ a £ 6 



- tu 2 + s 2 t 



3p(l) 
-*2 



X 



1357rM 3 s 2 t(s + m) 4 
| - u 2 (jt 2 + 12*u + 6m 2 ) - 12su(t + uf + s 2 (-7t 2 - 24tM - 



- 12s 3 (t + u) - 6s 4 



- u 2 (ht 2 + 12tu + 6m 2 ) - 12sm 2 (t + u) + sH (5t + 12m) 



+ 12s 3 (t + u) + 6s- 



da 
dt 



(gq 



cc 



lo(8) 



5gt (O 



lo(8) 



1152vrMs 2 t(s + M) 



x <^m 2 + s 2 + 



2 i 2 
M + S 



3^(8) 



5 



( ?s 6 (o[3 1 s{ 8) ])(4s 2 -sm + 4m 2 ) 
1728ttM 3 s 3 m(s + m) 2 



do" / 



x I - m 2 - 2tw - 2t 2 - 2st - s 2 + £ a £ 6 



- u (2t + u) + 2st + s 2 



lp(8) 



q 



967rM 3 s 2 t(s + m) z 



X <U 2 + S 2 + 



2 i 2 
— M + S 



da 
dt 



(gq -> cc 



3p (8) 







3p(8) 
-in 



(3s + 2t + 3m) 2 



864vrM 3 s 2 t(s + u) 



x {u 2 + s 2 + 



2 i 2 
— M + S 



da / 

dt ^ C " 



3p(8) 



-5^ 6 (g [ 3 A (8) 

" 432ttM 3 s 2 (s + m) 
x J tM 2 + 4sm (t + u) + s 2 (t + 4m) + f a f 6 



- tM 2 + S 2 t 



d<7 
dt 



(#3 



cc 



Sp(8) 



3p(8) 
-< 2 



4327rM 3 s 2 t(s + m) 
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X 



l - u 2 (jt 2 + 12tu + 6m 2 ) - 12su(t + uf + s 2 (-7t 2 - 2Uu - 12m 2 ) 



12s 3 (t + u) - 6s 4 + Ub 



u 



(ht 2 + \2tu + 6m 2 ) - 12sm 2 (t + u) 



+ s 2 t (5t + 12m) + 12s 3 (t + u) + 6s 4 



do" / 



io(i) 



= 



)(t 2 +M 2 ) 



277rMs 3 (t + M)^ 



da / _ 
dt 



cc 



3p(l) 



3 Po (1) 



(t 2 + m 2 ) (2s + 3t + 3m) 2 



8l7rM 3 s 3 (t + m) 



da 



6 I <n TspC 1 ) 



dt 



x 



gg — > cc 



Sp(D 



8^(0 ^ 



do" / 

-(gg-cc 



4tM (t + u) + s (t 2 + Atu + m 2 ) 

-8g° (O 



8l7rM 3 s 2 (t + m) 4 



3p(l) 
*2 



3 A (1) 



x 



405vrM 3 s 3 (t + uf 
6{t + m) 2 (t 2 + m 2 ) + 12s (t 3 + 2t 2 u + 2tM 2 + m 3 ) + s 2 (7t 2 + 12tu + 7m 2 



do" / 
-(gg-cc 



lo(8) 



= 



x -l+Ub>, 



da 
dt 



(gg - cc [ 3 Sf } ] <?) 



-5gt(Q [^ 8) ])(t 2 + M 2 ) 
4327rMs 3 (t + M) 2 



g 6 s (0 fe{ 8) l)(4t 2 -tM + 4M 2 ) 



6487rM 3 s 2 tM(t + uf 



x 



t 2 + u 2 + 2s (t + u) + 2s 2 ^ - 1 + >, 



da 



(gg 



cc 



lp(8) 



[l P f)|\ (t2+M2) 



367rM 3 s 3 (t + M)^ 



x -1 + Ub>, 
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da 
~dt 




5g« (O [ 3 P (8) ] ) (t 2 + u 2 ) (2s + 3t + 3u) 2 



(A.45) 



da 
~dt 



(A.46) 



da 
dt 



x [6(t + m) 2 (t 2 + u 2 ) + 12s (t 3 + 2t 2 u + 2tu 2 + u 3 ) + s 2 {jt 2 + I2tu + 7m 2 ) 



(A.47) 
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Figure 1: The unpolarized cross sections da/dpr and da/dy (first panel) and the double 
longitudinal-spin asymmetry All (second and third panels) of pp — > J /if) + X at RHIC- 
Spin, with y/S = 200 GeV, are studied as functions of pr (left column) and y (right 
column) in NRQCD and the CSM. The shaded bands indicate the theoretical uncertainties 
in the NRQCD and CSM predictions evaluated with the default PDFs. In the case of 
All, these uncertainties are compared with the spread due to a variation of the polarized 
PDFs. 
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Figure 2: Same as in Fig. ^ but for 7D — > J/ip + X in Experiment E161, with = 
45 GeV. Inelasticity z is used instead of rapidity y. In the first panel, the resolved-photon 
contributions are also shown separately. 
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Figure 3: Same as in Fig. ^ but for e + e — > e + e J/ij) + X in the e + e mode of TESLA, 
with v 7 ^ = 500 GeV, P(e+) = 60%, and P(e~) = 8( 



/o. 
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Figure 4: Same as in Fig. d but for e + e — > e + e J/ip + X in the 77 mode of TESLA, 
with VS = 500 GeV, P(e~) = 80%, and P(laser) = 100%. 
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Figure 5: The unpolarized cross section B(J/ip —>■ fi + {i )d 3 a/dy d 2 px of pp — ► J/ip + X 
followed by J/tj> -»■ //+/^ measured by PHENIX [22^ at RHIC, with v 7 ^ = 200 GeV, in 
the interval 1.2 < y < 2.2 as a function of is compared with the NRQCD and CSM 
predictions. 
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